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1 . INTRODUCTION 

Burbidge et al (1957) and Cameron (1957) laid out the framework for our 
understanding of the formation of the heavy nuclei (those nuclei with mass 
number A £ 70). From systematics in the solar system abundance distribution, 
Burbidge et al determined that the heavy nuclei were formed in three distinct 
nucleosynthetic processes, which they termed the r-, s-, and p-processes. That 
we still use these terms today is a credit to the soundness of this work done 37 
years ago. 

We may understand how Burbidge et al and Cameron arrived at their con- 
clusions from Figure 1. One population of nuclei, the s-nuclei, shows an 
abundance distribution with peaks near mass numbers 87, 138, and 208. These 
nuclei are made in a slow neutron-capture process, the s-process. A rapid 
neutron-capture process, the r-process, is responsible for the r-nuclei, whose 
abundance distribution shows peaks at mass numbers 80, 130, and 195. The 
p-process is responsible for production of the rarer, more proton-rich heavy 
isotopes (the p-nuclei) that cannot be made by neutron capture. 

The first quantitive evaluations of the ideas of Burbidge et al and Cameron 
came to light in the early 1 960s with work on the s-process (Clayton et al 1 96 1 , 
Seeger et al 1965) and the r-process (Seeger et al 1965). These calculations 
further elucidated the mechanisms for heavy-element formation and showed 
the plausibility of the framework developed in the 1950s. Subsequent work 
has focused on determining the astrophysical sites where the r-, s-, and p- 
processes occurred with the help of improved nuclear details, stellar models, and 
abundances. A goal of this paper is to review the recent progress astrophysicists, 
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Figure 1 The solar system abundances of r-nuclei, s-nuclei, and p-nuclei, relative to Si = 10 6 . 
Only isotopes for which 90% or more of the inferred production comes from a single process are 
shown. The data are from Anders & Grevesse (1989) and Kappeler et al (1989). 


astronomers, and physicists have made in these directions and to point out the 
problems that remain in our understanding of the formation of the heavy nuclei. 
Another, perhaps deeper, goal is to to seek some understanding of why there 
are three major processes available to nature for synthesis of heavy elements. 

It is impossible for a single paper to cover all relevant aspects of the r-, s-, and 
p-processes; therefore, where possible, references to other reviews are given. 
Readers should turn to these reviews for more details. Nevertheless, it is hoped 
that the present paper gives some flavor for the rich field of heavy-element 
synthesis. 

2. GENERAL CONSIDERATIONS 

The Master said, “Ssu, I believe you look upon me as one whose aim is 
simply to learn and retain in mind as many things as possible.” He replied, “That 
is what I thought. Is it not so?” The Master said, “No; I have one thread on which 
I string them all.” 


The Analects of Confucius 
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The r-, s-, and p-processes are distinct nucleosynthetic mechanisms. They occur 
in different environments and under quite different conditions. Nevertheless, 
it is useful to seek some unifying concept by which we may understand these 
disparate processes. We will see that entropy is the concept we need. Careful 
consideration of entropy in the various nucleosynthetic processes will clarify 
our discussion and give us insight into how these processes occur and why 
they occur where they do. With this insight we will see how the r-, s-, and 
p-processes are each unique answers to the same question: “How does nature 
produce heavy elements?” 

2.1 Entropy and Equilibrium 

Let us begin our discussion by considering a given thermally isolated system at 
constant volume. The system has total energy Eq. The entropy of this system 
is 


S = k InT, (1) 

where k is Boltzmann’s constant and V is the number of energetically allowed 
macroscopic states available to the system. By macroscopic state we mean 
a particular distribution of the constituents of the system among their single- 
particle quantum-mechanical states. In essence, a macroscopic state of the 
system is one particular way the constituents of the system can share the total 
energy Eq. Suppose the system does not have all macroscopic states of energy 
Eo available to it. In this case, the entropy is less than its maximum. The 
system will evolve by the Second Law of Thermodynamics and add more 
macroscopic states to its repertoire. In this evolution, the entropy will thus 
increase. The system will continue to evolve until all macroscopic states of 
energy Eo are available to the system. Once the system reaches this point, it 
is at maximum entropy and experiences no further evolution. The system has 
attained equilibrium. We thus see that equilibrium, or maximum entropy, is the 
evolutionary endpoint of any thermally isolated system. 

We can now use these considerations to ask what happens to the nucleons 
and nuclei in some nucleosynthetic environment. If the system of nucleons 
and nuclei is thermally isolated and out of equilibrium, it will evolve towards 
equilibrium. Given enough time, the system will reach equilibrium and attain 
maximum entropy. At this point, it is in nuclear statistical equilibrium (NSE), 
and it is a simple matter to compute the abundance of any nuclide (Burbidge 
et al 1957). We find (e.g. Meyer 1993) that, for a nuclear species of atomic 
number Z and mass number A, the abundance per baryon K(Z, A) is 

Y(Z, A) = G(Z , A)[<(3) /1 - , 7r (, - / ' )/2 2 <3/ '“ 5)/2 ]/4 3/2 (-^L) 3</ ' 


( 2 ) 
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where G(Z, A) is the nuclear partition function, f (3) is the Riemann zeta 
function of argument 3, T is the temperature, m n is the mass of a single baryon, 
(p is the photon-to-baryon ratio, Y p is the abundance per baryon of protons, K n 
is the abundance per baryon of neutrons, and B{Z, A) is the binding energy of 
nucleus (Z, A). We note that <p is given by 


<P = 


2 1 f(3 )(kT) 3 

n 2 (tic) 3 p'Na 


( 3 ) 


where is Avagadro’s number and p is the baryon mass density. The binding 
energy of nucleus (Z, A) is 


B(Z t A) = [Zm p + Nm n - m(Z , A)]c 2 , 


( 4 ) 


where N = A — Z and m(Z , A), m p , and m n are the masses of nucleus (Z, A), 
the proton, and the neutron, respectively. 

From Equation (2), we see that the NSE abundance of nuclei is nonzero. 
This might at first be surprising. If we combine free neutrons and protons into 
nuclei, we decrease the number of free particles of the system. This would 
yield fewer ways of sharing the total energy of the system and thus decrease the 
number of macroscopic states available to the system. The nuclear reactions 
also release binding energy, however, which increases the number of photons 
in the system, the energy available to leptons, and the excitation energy in the 
nuclei. These effects increase the number of ways the system can share the 
total energy of the system and, hence, can increase the number of macroscopic 
states available to the system. This increase can more than compensate for the 
decrease in the number of states due to the loss of free particles and can lead 
to an increase in the entropy. Once the system has evolved to the point that it 
experiences no net increase in the number of macroscopic states by changing 
the abundance of any particular nucleus, the system has reached NSE. 

Which nuclei dominate the abundance distribution in NSE? This depends 
on the photon-to-baryon ratio or, equivalently, the entropy per baryon since 
this latter quantity scales monotonically with the photon-to-baryon ratio (e.g. 
Meyer & Walsh 1993). From Equation (2) we see that the abundance of some 
heavy nucleus (Z, A) depends on (p [ ~ A . At some given temperature then, the 
larger <p is, the smaller will be the abundance of the heavy nucleus (Z, A). 
Fewer heavy nuclei means more light nuclei and free nucleons. The strong 
dependence of the NSE abundances on (p will be crucial for the r-process. This 
dependence on (p is apparent in Figure 2. If we fix the temperature, the larger 
(p is, the more likely nucleons are free or contained in light nuclei. 


2.2 Equilibrium Nucleosynthesis 

We have seen that the goal of any nucleosynthetic process is equilibrium. How 
does this goal vary with temperature? The answer is important because nature 
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Figure 2 The dominant species in NSE for a gas with electron-to-baryon ratio Y e = 0.5. The 
different regions of the plot show at what temperatures and photon-lo-baryon ratios 0 the various 
species dominate the gas. At high temperatures and high photon-to-baryon ratios, free neutrons and 
protons (labeled n, p) dominate. At lower temperatures and photon-to-baryon ratios, the nucleons 
are mostly locked up in 4 He nuclei. At still lower temperatures and photon-to-baryon ratios, the 
gas will predominantly be a distribution of iron-group nuclei. 


synthesizes nuclei at high temperature and ejects them into the medium between 
the stars at low temperature. If the nuclei stayed in equilibrium throughout this 
processing, we could use Equation (2) to determine the abundance of the various 
species. That the nuclei emerging from some nucleosynthetic process are not 
in general in an equilibrium distribution gives us important clues about how 
they were in fact made. We will turn to this point after discussing equilibrium 
nucleosynthesis in some detail. 

Consider a system in NSE at some temperature. If the system expands 
adiabatically, the temperature will decrease. By Le Chatelier’s Principle we 
know the system will respond to this stress by tending to counter it. How does 
the system do this? The answer is that nuclear reactions occur and release 
binding energy. This energy heats up the system and tends to counter the 
temperature decrease from the expansion. Via these nuclear reactions, the 
system finds a new equilibrium if the reactions occur faster than the expansion. 

As we have noted in Section 2.1, if nuclear reactions assemble free particles 
and light nuclei into heavier nuclei, there is a loss of macroscopic states and 
a decrease in the entropy in baryons (free nucleons and nuclei). The binding 
energy released in the reactions goes into photons and leptons. This increases 
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the numbers of macroscopic states and, hence, the entropy in the photons and 
leptons. We see that in building up the abundances of heavy nuclei, we transfer 
entropy from the baryons to the photons and leptons. This is apparent in 
Figure 3. Suppose a system at an entropy per baryon of 10£ begins in NSE at 
T 9 = 10(79 = 77 10 9 K) and expands adiabatically. The system cools. Because 
the system maintains equilibrium, entropy leaves the baryons as the baryons 
assemble themselves into heavier nuclei via nuclear reactions. This entropy 
goes into the photons and leptons (electron-positron pairs). 

The entropy in baryons continues to fall until T 9 ~ 4. At this point, the 
system has released essentially all its nuclear binding energy. There is little 
further transfer of entropy from the baryons to the other constituents of the 
system. We also note that beginning around 7g = 4 the electron-positron pairs 
annihilate into two photons. This transfers entropy from the — e~ pairs 
into the photons adiabatically. At late times, after the — e~ pairs have 
annihilated, the system is left with a residual abundance of electrons to ensure 
charge neutrality. After the pairs have disappeared, the system is a mix of 
relativistic photons and nonrelativistic nuclei and electrons. Now for adiabatic 
expansion, p oc T 3 for relativistic particles and p <x T 3/2 for nonrelativistic 
particles. The system at late times, because it is a mix of these two types of 
particles, actually expands such that p a T h where b is between 3/2 and 3. 
This means that entropy now is transferred from the photons to the translational 



Figure 3 The entropy per baryon in the consitituents of a system expanding adiabatically with 
total entropy per baryon 10* and electron-to-baryon ratio Y c . Note that the sum of the entropies of 
the constituents is always 10&. 
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entropy of the nonrelativistic particles. This is why the entropy in photons 
declines after T 9 = 1 and that in leptons and baryons rises. 

What nuclei are made in such an expansion? Again we note that nuclear 
reactions tend to occur only if there is a release of binding energy to compensate 
for the loss of macroscopic states due to the decrease in the number of free 
particles. As our system expands and cools, the nuclei present in equilibrium 
will increase in mass. They will continue to increase in mass until they have 
reached the nuclide with the largest binding energy per nucleon. Once this 
nuclide dominates the abundance of nuclei, there is no further evolution in the 
abundances. This is because if the system attempts to arrange the nucleons 
into a more massive nucleus with a lower binding energy per nucleon, there 
will be a net decrease in the total binding energy of the system. Reactions 
giving this rearrangement of nucleons lead to a decrease in the number of free 
particles and in the energy available to the photons and leptons. There will thus 
be fewer ways to share the energy of the system and a decrease in the entropy. 
An expanding and cooling system can maintain NSE by driving the system to a 
composition dominated by the nuclear species with the largest binding energy 
per nucleon. 

Figure 4 shows the mass fraction X(Z, A) — AY(Z , A) in NSE at tempera- 
tures 7g = 1-10 for an entropy per baryon of 1 0/c and a net electron-to-baryon 
ratio Y e = 0.5. If we imagine the system begins at T 9 = 10, alpha particles 



Figure 4 The mass fraction of species in NSE as a function of temperature for a system with 
entropy per baryon I0& and Y c =0.5. 
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dominate the mass of the system. As the system expands and cools adiabat- 
ically, the abundances of heavier nuclei increase. At late times, however, all 
nucleons are locked up in 56 Ni nuclei. This is because 56 Ni has one proton 
for every two nucleons ( Y c = 0.5), and it is the nucleus at Y c = 0.5 with the 
strongest binding energy per nucleon. 

Figure 5 shows what happens if the entropy per baryon is \00k instead of 
10 k. Because of the larger entropy, the system has more free particles and light 
nuclei at a given temperature. This is why free neutrons and protons dominate 
the abundances around T 9 = 10 in Figure 5 while alpha particles dominate the 
abundances in Figure 4. We also see that alpha particle abundances build up 
only after T 9 drops below 10 for entropy per baryon of 100&. As in the lower 
entropy case, however, the nucleons eventually end up all in 56 Ni. 

What happens now if Y c is different from 0.5? Figure 6 shows what happens 
if the entropy per baryon is 10 k and K e = 0.4. Here there are six neutrons for 
every four protons. We see that the system at late times locks up all its nucleons 
into 70 Ni. For K e = 0.4, this is the nucleus with the largest binding energy 
per nucleon. 

From Figure 6 we see that the final nucleus on which the system converges 
at low temperature is dependent on F e because the final nuclei present must 
accommodate all of the neutrons and protons initially in the system. This final 
nucleus, however, is the nucleus with the largest binding energy for that T e . 



Figure 5 The mass fraction of species in NSE as a function of temperature for a system with 
entropy per baryon 100 k and Y<. = 0.5. 
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Figure 6 The mass fraction of species in NSE as a function of temperature for a system with 
entropy per baryon 10& and Y e = 0.4. 


Figure 7 gives an idea of what that nucleus is for K e ranging from 0.35 to 0.6. 
In all cases, the nuclei are iron-group nuclei (Z = 26-34). 


2.3 How to Make Heavy Nuclei 

We are now in a position to ask the question “How does nature make heavy 
nuclei?” We have seen in the previous sections that if we eject nuclei into 
cold interstellar space, and if these nuclei are always in NSE, then these nuclei 
must be those with the largest binding energy per nucleon for whatever Y c is 
appropriate for the environment in which the nuclei find themselves. From 
Figure 7, however, we see that for a large range of Y e , these nuclei are simply 
iron-group nuclei. If NSE always pertains, stars and supernovae can only eject 
iron-group nuclei. Our observations of uranium on Earth tell us that this is not 
what happens. 

The escape from this dilemma is the fact that nucleosynthetic systems cannot 
always be in NSE. There are only two ways this can happen. The first possi- 
bility is that the system never has time to come into NSE before the star ejects 
the nucleons into the interstellar medium. In this case the nucleons assemble 
themselves part of the way up to iron. This is the “falling short of equilib- 
rium” scenario. The second possibility is that the system begins in NSE at high 
temperature. As the system expands and cools, the equilibrium changes. As 
the temperature drops, some nuclear reactions slow down. Eventually these 
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Figure 7 The proton number Z {top) and the mass number A {bottom) of the species with the 
largest binding energy per nucleon in each bin in Y c . Where experimental data are available, they 
are from Wapstra et al (1988). Where experimental data are unavailable, nuclear masses are from 
Moller& Nix (1988). 


reactions become too slow to allow the system to maintain equilibrium. This 
is the “freezeout from equilibrium” scenario. These are the only two ways a 
nucleosynthetic system can be out of equilibrium, and they are the only two 
options nature has to assemble heavy nuclei. 

How in fact does nature make heavy elements in these two scenarios? In the 
freezeout from equilibrium, a nucleosynthetic system falls out of equilibrium. 
At this point the composition reflects the last NSE abundance distribution the 
system attained before it fell out of equilibrium. This composition is a mix of 
free nucleons, light nuclei, and iron-group nuclei. While some of the reactions 
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necessary to maintain equilibrium become too slow, others such as capture 
of free nucleons and light nuclei on iron-group nuclei can continue. The iron- 
group nuclei can serve as “seeds” for the capture of the remaining free nucleons 
and light nuclei. These captures can produce heavy nuclei. Because the system 
assembles its own seed nuclei in such a scenario, the system can make heavy 
elements without any pre-existing seed nuclei. A process that does not require 
pre-existing seed nuclei to make new nuclei is called primary. 

In the other scenario — the falling short of equilibrium scenario — the material 
never achieves NSE because the timescale to reach equilibrium is always too 
long compared to the dynamical timescale of the system. The system assembles 
light nuclei into somewhat heavier nuclei but these are still less massive than 
the iron-group nuclei characteristic of NSE. An accidental effect of the main 
nuclear reactions between the light nuclei in this scenario is the liberation of 
nucleons. If iron-group or heavier nuclei already exist in the system, they can 
capture these nucleons to produce new heavy nuclei. In this scenario, then, the 
system must have pre-existing seed nuclei in order to produce heavy elements. 
We call a process that requires pre-existing seed nuclei to make new elements 
a secondary process. 

These are the only two ways nature can assemble heavy elements. We should 
not be surprised then that there are two major distributions of heavy nuclei — the 
r-process and the s-process distributions. We shall see that the r-nuclei in our 
Solar System likely formed in an environment that experienced a freezeout from 
equilibrium while the s-nuclei must have formed in an environment that was 
striving for, but never reached, NSE. The differing character of these scenarios 
results in the different character of the r- and s-process abundance distributions. 

Once an abundance of heavy elements is available, nature may make modi- 
fications to it by exposing it to a flux of photons, neutrinos, or nucleons. Such 
events are probably responsible for the production of the majority of p-nuclei. 

3. THE r- PROCESS 

Nature does not endure sudden mutations without great violence. 

Francois Rabelais, Gargantua 

We turn now to a discussion of the r-process. We seek to understand in some 
detail how and where the r-process occurs. How the r-process occurs depends 
on whether it is a primary or secondary process. We will see that the evidence 
available to us today indicates that the r-process is a primary, freezeout from 
equilibrium process. As for where the r-process occurs, the rapid timescales 
associated with the r-process point to violent events such as supernovae or 
disruptions of neutron stars. Winds from nascent neutron stars are currently 
believed to be the most probable site. Our discussion will be brief. For more 
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details, the reader should turn to the many excellent reviews in the literature 
(e.g. Hillebrandt 1978, Mathews & Cowan 1990, Cowan et al 1991). 

3. 1 The Primary r- Process 

Let us return to the notion of a freezeout from equilibrium. As a system in NSE 
expands and cools, the abundances shift to maintain NSE. We have seen that 
in doing so the entropy of the system moves from the baryons into the photons 
and leptons and the abundance of heavier nuclei grows at the expense of free 
nucleons and light nuclei. Eventually the temperature of the system is too low 
or the abundance of reactants too small for certain reactions to go fast enough to 
maintain NSE. These reactions freeze out. The first reactions to freeze out are 
charged-particle reactions. Neutron-capture reactions can continue, however, 
because they are not impeded by a nuclear Coulomb barrier. The nuclei present 
at the time of the freezeout of the charged-particle reactions then eventually 
capture the remaining neutrons. 

It we suppose that the system is quite neutron rich, many free neutrons 
should exist after freezeout of the charged-particle reactions. One might imag- 
ine that as we increase the neutron richness of the system, we simply have 
more neutron-rich isotopes of the nuclei present in the system. However, there 
is a limit to how neutron rich the nuclei can get. This is because the nuclei 
eventually encounter neutron drip, at which point the binding energy of the 
next neutron captured is negative. Once the nuclei reach neutron drip, they 
cannot contain any more neutrons. Thus, for sufficiently neutron-rich mate- 
rial, there may be many free neutrons at the time of charged-particle reaction 
freezeout. The system establishes an equilibrium between the neutron-capture 
(«, y) and neutron-disintegration (y, n) reactions. Beta decays then occur 
which increase the charge on the nucleus and allow further neutron capture. 
This phase of the freezeout from equilibrium in which only neutron-capture, 
neutron-disintegration, and beta-decay reactions occur is the r-process. 

How neutron rich must material be to undergo an r-process? We know that 
the r-process produces uranium (A = 238) from seed nuclei (A % 50-100); 
therefore, there must be roughly 100 free neutrons per seed nucleus at the 
time of charged-particle reaction freezeout. The required neutron richness thus 
depends on the abundance of seed nuclei which in turn depends on the entropy 
per baryon in the system. 

If the entropy per baryon is less than or roughly 10&, NSE favors iron- 
group nuclei. These are the seed nuclei for the r-process once the charged- 
particle reactions have frozen out. Because the r-process needs around 100 free 
neutrons per seed, a seed nucleus like 78 Ni, typical in a low entropy, neutron-rich 
freezeout, requires a total neutron-to-proton (n/p) ratio of (50 + 100)/28 ^5.4 
or a Y e = p/(n -bp) ^ 0.16. This is quite neutron-rich material. 

If the entropy per baryon is high (> 100A:), NSE favors 4 He at high tern- 
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perature. As the temperature drops in the expanding material, NSE begins to 
favor iron-group nuclei. For large entropy per baryon, this occurs late in the 
expansion. At this late time, the two reaction sequences that begin the assembly 
of alpha particles into iron-group nuclei, 4 He + 4 He + 4 He — ► 12 C and 4 He + 
4 He + n -> 9 Be followed by 9 Be + 4 He -* 12 C + n, may not be operating at 
a significant level or may even have frozen out. We note that these reaction 
sequences rely on three-body interactions which are highly sensitive to the den- 
sity. The higher the entropy per baryon at a given temperature, the larger the 
photon-to-baryon ratio and the lower the density. A lower density results in 
slower three-body reaction rates. Thus the higher the entropy per baryon, the 
higher the temperature at which these three-body reactions freeze out and the 
lower the abundance of seed nuclei. Unlike the three-body reactions, alpha cap- 
tures on ,2 C and heavier nuclei occur rapidly and build up heavy nuclei. As the 
system falls out of equilibrium, there will be more alpha particles around than 
there would be in NSE. Some of these alpha particles capture on the iron-group 
nuclei present to make heavier nuclei (A ~ 100). Eventually these reactions 
freeze out also, and the system is left with an abundance of free neutrons, seed 
nuclei, and many 4 He nuclei. At this point the free neutrons can capture on the 
seed nuclei, but not on the 4 He nuclei. In this way an r-process occurs. 

The degree of neutron richness necessary for an r-process in a high entropy 
environment is less than in a low entropy environment. A typical composition 
for an entropy per baryon of 400A: and a neutron-to-proton ratio of 1.6 (K c = 
0.385) is 20% free neutrons by mass, 10% seed nuclei, and 70% alpha particles 
(Woosley & Hoffman 1992). Suppose the seed nucleus is Z = 35, A = 100 
(also typical), then the free neutron-seed nucleus ratio is 200 in the case: more 
than sufficient for an r-process. We thus see that the degree of neutron richness 
we need for a primary r-process directly depends on the entropy per baryon. 

We have to this point discussed how a primary r-process would occur in 
the freezeout from NSE in low entropy and high entropy environments. The 
r-process might also occur in a very low entropy or zero-entropy environment. 
Matter inside a neutron star is extremely neutron rich and highly-degenerate. 
If the neutron star is more than a few hours old, it is cold (T ^ 10 9 K) on a 
nuclear energy scale (e.g. Baym & Pethick 1979), implying that the entropy 
per baryon is quite low (;S 0.5k). At densities below nuclear density (p^2x 
10 14 gcm -3 ) but above neutron-drip density (p^4x 10 u gcm~ 3 ), extremely 
neutron-rich nuclei exist in strong equilibrium with degenerate neutrons and 
weak equilibrium with the degenerate electrons (Baym et al 1971). Above 
nuclear density, the material is comprised of free nucleons and electrons. Weak 
equilibrium forces the Y e of the material to be of order 0.05 or less (Lattimer et 
al 1985, Lattimer & Swesty 1991). 

Only the strong gravity of the neutron star keeps such matter from exploding 
apart. If a piece of cold neutron-star matter were to escape from the neutron star 
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into interstellar space, it would decompress. If this escape occurred without 
too much violence, there would not be a dramatic increase in the entropy. The 
material would remain cold. The material would consist of neutrons and nuclei, 
and as the material expanded, an r-process could occur (Lattimer et al 1977, 
Meyer 1989). In such a system, the only heating that would occur would be from 
the beta decays and nuclear fissions during the r-process, which are irreversible 
processes. Notice that such an r-process would not be a freezeout from NSE. 
The system was not in nuclear statistical equilibrium prior to expansion, and it 
would not attain NSE later in the expansion unless the beta decays and nuclear 
fissions that occur could drive the temperature up high enough, something 
that probably does not happen (Meyer 1989). On the other hand, this would 
be a freezeout from weak equilibrium. Thus, such an r-process would be 
primary because it is the formation of the neutron star that creates the seed 
nuclei and the weak equilibrium would erase the entire previous history of 
the nucleons. Moreover, material that began at densities above nuclear matter 
density would experience a phase transition from free nucleons into neutrons 
and nuclei during the expansion. In this case, the seed nuclei would form during 
the decompression. 

Now that we understand how the r-process occurs in freezeout from equi- 
librium, we may consider the question of what astrophysical sites could give 
low, very low, or high entropy r-processes. Let us begin with low entropy sites. 
The earliest site considered for a low entropy r-process was at the mass cut of 
a type II supernova (Burbidge et al 1957, Cameron 1957). The mass cut is the 
boundary between the matter that escapes into space from the supernova and 
the matter that remains as part of the remnant neutron star. The first serious 
time-dependent calculations of such an r-process were made by Seeger et al 
(1965), although this work assumed constant temperature and neutron number 
density. Later calculations treated the r-process as dynamical, that is, with vary- 
ing temperature and neutron density (Cameron et al 1970, Schramm 1973, Sato 
1974, Kodama & Takahashi 1975, Hillebrandt et al 1976, Hillebrandt 1978). 
While some of these workers found fairly good fits to the solar system r-process 
distribution, the models studied provided no natural way of explaining why the 
particular r-process distribution we see should emerge from a supernova and, 
more seriously, why a supernova should eject only a small mass of r-process 
matter (see Section 3.3). It is also important to note that current type II super- 
nova models do not yield K e as low as 0. 1-0.2 at the mass cut, as required by a 
low entropy r-process (e.g. Wilson & Mayle 1993, Woosley et al 1994). 

Due to the difficulties with the mass-cut site, astrophysicists turned to other 
supernova scenarios. In particular, workers sought other means of ejecting low 
entropy, neutron-rich r-process matter. Rotating stellar cores with (LeBlanc & 
Wilson 1970, Meier et al 1976, Muller & Hillebrandt 1979) and without mag- 
netic fields (Symbalisty 1984, Symbalisty et al 1985) can eject some neutron- 
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matter. It is not clear, however, that the supernova cores can attain the high 
rotation rates and magnetic fields required for such ejection. Moreover, there 
is no natural explanation for why the material has just the right conditions to 
make a solar system r-process distribution. 

What about very low entropy r-processes? Lattimer & Schramm (1974, 
1976) considered the tidal disruption of a neutron star by a black hole. The 
ejected neutron-star matter could then undergo a very low entropy r-process 
(Lattimer et al 1977, Meyer 1989). Neutron star-neutron star collisions could 
also lead to the ejection of neutron-rich matter (Symbalisty & Schramm 1982, 
Eichleretal 1989; see also Kochanek 1992 and Colpi etal 1989, 1991, 1993). If 
this material were not strongly disturbed during the collision, it would undergo 
a very low entropy r-process. It is likely in such an event, however, that the 
material would be shocked to entropies of order several k per baryon, so that 
in fact a low entropy would probably ensue (Evans & Mathews 1988). In any 
case, these sites suffer the two difficulties of uncertain occurrence rates and the 
lack of any natural reason a solar system r-process distribution should result. 

This leaves the high entropy r-process. A promising site for such an r-process 
is in the neutrino-driven winds from nascent neutron stars. We consider this 
site in more detail in Section 3.4. 

3.2 The Secondary r~Process 

A system that does not achieve equilibrium may also produce r-process ele- 
ments. As we have seen, as the bulk of the nuclei strive to reach the iron 
group, the reactions that carry them in this direction may release neutrons that 
then capture on pre-existing seed nuclei. If the number of neutrons released is 
sufficiently large, a solar system r-process abundance distribution may result. 
The nuclear reactions that would be the dominant producers of neutrons are 
(a, n) reactions, and in particular the reactions ,3 C(a, n) 16 0, 22 Ne(», «) 25 Mg, 
and 25 Mg(a, rc) 28 Si. 

A crucial feature of a secondary r-process is that it does not achieve («, y) — 
(y, n) equilibrium. Such equilibrium occurs when the flows from (n, y) and 
(y , n) reactions come into balance. The most abundant isotope of some element 
in such equilibrium is then the one for which the rate for that isotope (Z, A) 
to capture a neutron is equal to the rate for the resultant isotope (Z, A -j- 1) to 
suffer a disintegration (y, n) reaction. The neutron-separation energy S n , that 
is, the binding energy of the least tightly bound neutron, of this nucleus is, upon 
neglect of the A dependence of nuclear partition functions (e.g. Sato 1974, see 
also Howard et al 1993, and Meyer 1994), 

S n (MeV) = -^-(34.08+ 1.5 log 10 r 9 - 1.5log 10 « n ), (5) 

where n n is the neutron number density in units of cm -3 . In secondary r-process 
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sites, the neutron number density is typically ~ 10 19 cm~ 3 and 7^ ~ 1 ; thus, the 
dominant nuclei in (n, y) — (y,n) equilibrium have S n ~ 3.0MeV. Figure 8 
shows the neutron-separation energies for the isotopes of neodymium (Z = 60). 
We see that the dominant isotope in («, y) — (y, n) equilibrium should be the 
one with neutron number N & 110. The most neutron-rich beta-stable isotope 
of neodymium is l50 Nd (N = 90). The neutron sources in a secondary r-process 
would have to supply and maintain some 20 neutrons per seed nucleus in order 
to establish (n, y) — ( y , n ) equilibrium. This is something they cannot do. 

Let us now consider possible secondary r-process sites. One possible site is 
the helium shell of an exploding massive star (Truran et al 1978, Thielemann 
et al 1979; see also Cowan et al 1980, 1983, 1985; Blake et al 1981; Klapdor 
et al 1981). As the supernova shock wave traverses this shell, it heats up the 
material to a temperature Tg & 1 . Although this temperature is not high enough 
to force the material into NSE, it is high enough to drive material strongly in 
that direction. Among the main nuclear reactions that occur are (a, n) reactions 
that liberate neutrons which can then capture on the pre-existing seed nuclei. 
Truran et al (1978) found that neutron captures could modify an s-process seed 
abundance distribution into an r-process distribution. Later work indicated 
that a seed distribution that is enhanced with respect to the solar system heavy- 
element distribution is required to produce the solar r-nuclei. Such an enhanced 
distribution could result from s-processing prior to shock passage if protons 
were to mix down into the helium shell to make l3 C from the abundant 12 C [via 
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Figure 8 Neutron-separation energies for the isotopes of neodymium. The data are derived from 
Moller & Nix (1988). 




r-, s-, AND p-PROCESSES 


169 


12 C(/?, y) 13 N(/? + ) ,3 C]. It is apparent today, however, that the amount of 13 C 
required is unrealistically large (Cowan et al 1985, Cameron et al 1985), 

Lee et al (1979) considered supernova shock passage through the carbon 
shell. Here the 22 Ne(a, y) 25 Mg reaction produces the neutrons. 25 Mg is also a 
neutron poison, however; it absorbs many of the liberated neutrons before they 
have a chance to capture on heavier seed nuclei. Wefel et al (1981) found that 
some heavy nuclei could be synthesized, but not the bulk of the r-nuclei. 

Another intriguing secondary site is again the helium burning shell, but now 
the effects of neutrino inelastic scattering on 4 He nuclei are included. The 
neutrinos come from the cooling nascent neutron star resulting from the core- 
collapse event. These neutrinos spall neutrons from the 4 He (and other light) 
nuclei. The neutrons can then capture on the seed nuclei and drive an r-process 
(Epstein et al 1988). A detailed study of this “u-process” showed that some 
interesting nucleosynthesis may occur in such an event, but it could not have 
been a major contributor to the solar system r-process abundances (Woosley et 
al 1990). 

All of the secondary r-process models studied to date have had profound 
difficulties which have rendered them implausible. Moreover, as shown in the 
next section, there are other reasons for favoring a primary over a secondary 
r-process. Nevertheless, the study of secondary r-process models has been 
valuable because they may have important applications to isotopic anomalies 
in meteorites (e.g. Clayton 1989, Howard et al 1992, Cameron et al 1993). 

3.3 Observational Constraints 

We have discussed how the r-process works in both primary and secondary 
scenarios and possible sites for the r-process. Let us turn now to the “observa- 
tional” data to evaluate the plausibility of the proposed sites. These data will 
be astronomical, meteoritical, and nuclear. 

We begin with the observations of old stars carried out by Sneden and col- 
laborators (Sneden & Parthasarathy 1983, Sneden & Pilachowski 1985, Gilroy 
et al 1988; see also Wheeler et al 1989). These workers find that the elemental 
abundances, i.e. abundances as a function of Z, the proton number, of heavy 
elements in the atmospheres of these old stars match rather well the solar sys- 
tem r-process elemental abundances; s-process elemental abundances do not 
provide a good fit. One conclusion we may draw from this is that the r-process 
dominated the s-process at early times in our Galaxy. This would strongly sug- 
gest that the r-process is primary because a secondary r-process must produce 
r-nuclei from an already abundant population of s-nuclei. Another important 
conclusion is that the r-process mechanism has remained essentially the same 
throughout the Galaxy's history. A caveat to both of these conclusions is that 
they are based on elemental abundances, whereas a truly definitive identifica- 
tion of the heavy nuclei in old stars as r-process in origin would require isotopic 
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abundances, information that spectroscopy cannot yet generally give us. Nev- 
ertheless, the evidence for a primary r-process that has operated in the same 
fashion throughout the Galaxy’s history is compelling. 

The next item of evidence we may consider is the question of the timescale for 
heavy element formation. Mathews & Cowan (1990) and Mathews et al (1992, 
1993) have fit chemical evolution models to observations of the elemental abun- 
dances of europium (a mostly r-process element) and iron in the atmospheres of 
old stars. The only models that give good fits to the data are a primary r-process 
in type II and lb (that is, core-collapse of a massive star) supernovae or a sec- 
ondary r-process in which the neutron source itself is primary. An example of 
a primary neutron source in a secondary r-process event is the 13 C produced 
by mixing of protons down into the helium shell for supernova shock-induced 
helium or carbon burning. This is effectively a primary source because the star 
itself constructs the 13 C from its initial supply of protons and from 12 C made 
during helium burning. Another example is the neutrino- induced r-process. 
The neutrons in this source come mostly from inelastic scatterings of neutrinos 
on 4 He nuclei in the helium shell. The star itself produced almost all of the 
4 He nuclei, so the neutron source is primary. Other proposed sites seemingly 
ruled out by these timescale arguments are: tidal disruptions of neutron stars 
by black holes, neutron star-neutron star collisions, a secondary r-process with 
a secondary neutron source such as 22 Ne, neutron-star accretion disks (Hogan 
& Applegate 1987), core helium flash in low-mass stars (Cowan et al 1982), 
and classical novae (Hoyle & Clayton 1974). 

Experimental nuclear physicists have also provided us with important clues 
about the r-process. In particular, K.-L. Kratz and collaborators have measured 
the beta-decay lifetimes of r-process “waiting pointing” nuclei on the N — 
A - Z = 50 and N = 82 closed neutron shells (Kratz et al 1988, 1990). 
These nuclei have particularly strongly bound valence neutrons. Now the total 
abundance of the element Z is a sum over the abundances of all of the isotopes 
of that element: 


y z = J2 y ( z . Ay (6) 

A 

The average rate of beta-decay flow out of element Z is 

Xz = (7) 
Y z 

where X^(Z, A) is the beta-decay rate of nucleus (Z, A). The abundance Y z 
decreases during the r-process by beta decay out of Z into Z -f 1 and increases 
by beta decay into Z from Z — I ; thus, 


dY z 

dt 


zYz + ^z-i Yz- 1. 


( 8 ) 
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Where the r-process path crosses the N = 50 and N = 82 closed shells, a 
single isotope dominates the abundances because of the strong neutron binding. 
For these nuclei, Y z ~ Y(Z, Z + 50) for the N = 50 closed neutron shell 
and Y z & T(Z, Z + 82) for N = 82. If the conditions for («, y) — (y, n) 
equilibrium hold for a long enough time, the system will achieve approximate 
steady beta flow such that dY z /dt — > 0. In this case k z Y z — ► constant. 
What is particularly significant is that the decay rates found by Kratz and 
co-workers times the appropriate abundances along the closed neutron shells 
as inferred directly from the solar system r-process data show that A Z Y Z is 
approximately constant for the N — 50 and /V = 82 waiting-point nuclei. 
This is strong, albeit circumstantial, evidence that the r-process achieved the 
conditions of (n, y) — (y, n) equilibrium and steady beta flow. This would 
indicate that the r-process is primary since a secondary r-process does not 
achieve (n,y) — (y, n ) equilibrium. Another important result found from this 
work is that the constant k z Y z for the N = 50 closed neutron shell is different 
from that for the N = 82 closed shell. The r-process did not achieve global 
steady beta flow (Kratz et al 1993). We interpret this to mean that the r-process 
abundance distribution actually results from a sum of components resulting 
from a varying set of r-process conditions — no single set of conditions gives 
rise to the solar system abundance curve. While this point was long appreciated 
by theoreticians (Seeger et al 1965, Kodama & Takahashi 1975, Hillebrandt et 
al 1976), the work of Kratz and collaborators establishes this result directly 
from nuclear physics and meteoritic abundances. The actual r-process site in 
nature must naturally give rise to a varying set of r-process conditions. 

The final observational constraint we consider is the mass of r-process mate- 
rial in the Galaxy. From meteorites and abundances in the Sun we can infer that 
the mass fraction of r-nuclei in the Galaxy is 2 x 10 -7 (from the data in Anders 
& Grevesse 1 989). If the mass of our Galaxy is 1 .5 x 10 1 1 M 0 , there are some 
10 4 M g of r-process material in our Galaxy. Now the rate of supernovae in our 
Galaxy is between 0. 1 yr _1 and 0.01 yr -1 (e.g. Tammann 1982, Van den Bergh 
& Tammann 1991) and the Galaxy is of order 10 10 years old; therefore, there 
have been some 10 8 to 10 9 supernovae in our Galaxy’s history. If each of these 
supernovae produced r-process material, we expect each supernova to make 
10~ 5 to 10" 4 Af Q of r-nuclei. This is a tiny fraction of the total mass ejected in 
a supernova and its smallness provides an important constraint on the site of 
the r-process. 

From the astronomical and nuclear observations, we have found that the r- 
process is likely a primary event occurring in core-collapse supernovae. If the 
r-process comes from low entropy, neutron-rich material ejected from the core 
in a supernova, then we know from the discussion in Section 2.1 that we need 
material with K c ;$ 0.2 to make actinide nuclei. Such material must come from 
fairly deep in the core. Hillebrandt et al (1976) found that a supernova needs 
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to eject of order 0.1 A/q of material to get T e ^ 0. 1-0.2 material out of the 
core. Such a large amount would overproduce the r-nuclei by a factor of 10 3 if 
all supernovae ejected this much neutron-rich matter. It may be that only rare 
supernovae, such as those with high magnetic fields and rotation rates (LeBlanc 
& Wilson 1970, Meier et al 1976, Muller & Hillebrandt 1979, Symbalisty 
1984, Symbalisty et al 1985), can eject neutron-rich matter. It remains to be 
seen, however, whether the high magnetic fields and/or rotation rates in these 
scenarios are indeed achieved in nature. In conclusion, we have no natural 
scenario for a low entropy primary r-process that yields 10~ 4 A/ o of r-process 
material per event. Fortunately, we do have a high entropy r-process that gives 
this yield of r-nuclei per supernova. We turn to this scenario in the next section. 

3.4 The r-Process in Nascent Neutron Star Winds 

Woosley & Hoffman (1992) proposed the winds from nascent neutron stars as 
the site of the r-process. In this section we discuss the physics of these winds. 
Only detailed models will allow us to determine whether these winds are indeed 
the site of the r-process. On the other hand, simple arguments will illustrate the 
basic features of these winds and why they are attractive as an r-process site. 

A core-collapse (i.e. type II or lb) supernova may leave a hot (kT & lOMeV) 
neutron star as a remnant (e.g. Bethe 1990). This neutron star cools by neu- 
trino emission on a timescale of order 10 seconds. Salpeter & Shapiro (1981) 
were among the first to consider the thermal evolution of such a young neu- 
tron star. They were able to show that, although the neutrino luminosity is 
always sub-Eddington, the photon luminosity can be super-Eddington. Thus, 
for sufficiently hot nascent neutron stars, there is a neutrino-driven wind. 

Duncan et al ( 1986) repeated the arguments of Salpeter & Shapiro and showed 
that a spherically symmetric, nascent neutron star of radius 10.6 km and mass 
1 AMq has a super-Eddington photon luminosity if the neutron star temperature 
is £ 0.4 MeV. Duncan et al also showed that the star cannot stabilize itself 
by changing its size or by transporting energy by convection. A wind must 
therefore blow from the surface of the neutron star. Moreover, Duncan et al 
showed that the mass loss rate in this wind is of order 10~ 5 M© s _1 , which gives, 
if the wind lasts for the early cooling time of the nascent neutron star (^10s), 
about 10 ” 4 Mq. If this material forms r-process nuclei, the wind naturally gives 
the correct amount of r-process matter per supernova. [For more details on this 
wind, see Woosley et al (1994) and references therein.] 

What about the entropy in the wind? The net heating of a matter element 
lifting off the surface of the neutron star is governed by the heating due to 
neutrino interactions with the wind material and cooling due to emission by the 
matter. The heating goes as F v o <£ v ) , where a is the neutrino-matter interaction 
cross section, F v is the number flux of neutrinos, and (E v ) is the average neutrino 
energy. F v oc L V (E V )~ { r~ l where L v is the neutrino luminosity and r is the 
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radial position of the matter element, a a (£ y ) 2 (e.g. Tubbs & Schramm 1975). 
The heating thus goes as L v (E v ) 2 r~ 2 . {E v ) and L v fall off on the neutron-star 
cooling timescale of roughly 10 seconds. The matter elements move out on 
timescales faster than 10 seconds, so the heating rate falls off roughly as r ~ 2 . 
The cooling goes as 7^, where T m is the local matter temperature. Above the 
surface of the neutron star, 7^ falls off more steeply than r -2 . As the mass 
element lifts off the star, the initial heating is slow because the matter and 
neutrino temperatures are nearly equal. As the mass element moves out, the 
matter temperature T m drops. Once the mass element passes the “gain radius,” 
where the heating and cooling rates are equal, T m is too low for the material to 
cool off as fast as it is heated (Bethe & Wilson 1985). As heat is added, the 
entropy rises. In this way, the entropy can reach values of 100/: or more before 
neutrino interactions with the wind material freeze out. These are the values 
required for a high entropy r-process. 

Note that the entropy in mass elements leaving the neutron star at late times 
will be larger than in mass elements leaving the star at early times. This is 
because most of the heating occurs fairly near the surface of the neutron star 
because the neutrino flux and hence the net heating rate falls off as 1/r 2 . As a 
neutron star ages, it shrinks in radius [from ~ 1 00 km at a few tenths of a second 
after core bounce to ~10 km at several seconds after bounce in the models of 
Wilson & Mayle (1993)]. The decrease in the initial r from which the mass ele- 
ments begin increases the heating rate more than the slow fall off in L v and (E v ) 
decrease it. The net heating and entropy in the later mass elements is thus larger. 

The last question we must consider is the neutron richness of the wind mate- 
rial. Y e is set in the wind by the reactions v e + « — ► p + e~ and v c + p — ► n + e^. 
If the fluxes and energies of the v e and u e were equal, Y c would be slightly larger 
than 0.5 because the mass of the proton is slightly less than that of the neutron. 
As a nascent neutron star cools, however, it becomes neutron rich. The opacity 
for v c s becomes larger than that for v e s because of the reaction v e +n — ► p + e ~ . 
The v e s thus have a longer mean free path and originate deeper in the neutron 
star. This means that they are more energetic than the u c s. This necessarily 
drives the material neutron rich. At late times (/ ~ 10 s after core bounce) 
(Ef, c ) ~ 2(E Vc ) which makes Y c ~ 0.33 (Qian et al 1993). This is certainly 
neutron rich enough for a high entropy r-process. 

As a last point, let us re-emphasize that the entropy and T e in a mass element 
vary according to when that mass element lifts off the neutron star. Each 
mass element thus undergoes a somewhat different nucleosynthesis. The final 
r-process abundance distribution, however, is a sum of all of these different 
components. The r-process in nascent neutron star winds thus naturally satisfies 
the requirement imposed by the work of Kratz et al (1993). Notice that the wind 
dynamics and thermodynamics are completely determined by the mass, radius, 
and temperature of the neutron star. Any given neutron star probably passes 
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through roughly the same sequence of temperatures and radii as a function of 
time as it cools. We thus expect to get essentially the same r-process out of 
every supernova. This satisfies our expectations from the observations of the 
r-process elements in old stars (see Section 3.3). 

Does a solar system distribution naturally emerge in such a wind? Meyer 
et al (1992), using a schematic model based on output from Wilson & Mayle 
(1993), found the resulting abundances matched the solar system distribution 
quite well (see also Howard et al 1993 and Takahashi et al 1994). A more 
detailed model, using mass element trajectories calculated directly in Wilson 
and Mayle’s supernova code produced abundances that also agree well with 
the solar distribution (Woosley et al 1994). This latter model also gives the 
correct amount of r-process mass per supernova (~ 10 _4 Af G ). Confirmation of 
nascent neutron star winds as the site for the r-process will require a full survey 
of the nucleosynthesis in detailed, realistic wind models. Nevertheless, nascent 
neutron star winds seem extremely promising as the site for the r-process. 


4. THE s-PROCESS 

Must we, as Solon advises, always keep the goal in sight? 

Aristotle, Nichomachean Ethics 

The s-process is the other major nucleosynthetic process that assembles heavy 
elements. We know that the s-process path in the neutron number-proton 
number plane crosses the neutron closed shells at the valley of beta stability. 
This tells us that the s-process occurred in an environment with a much lower 
neutron density than the r-process. Also, the s-process occurred over a much 
longer time period. 

In this section we seek to understand how the s-process occurs. We then turn 
to the question of s-process sites. Finally we consider constraints on those sites. 

4. 1 The s-Process Mechanism 

Because of the neutron densities and timescales inferred for the s-process from 
the abundance peaks, we can infer that the s-process is not a freeze out from 
equilibrium. Instead, it is a neutron-capture process that occurs in a system 
striving to reach equilibrium, but falling short of its goal. The main reactions 
carrying the bulk of the nuclei towards the iron group can liberate neutrons. 
Pre-existing seed nuclei capture these neutrons and produce the s-nuclei. The 
s-process is clearly a secondary process. 

The dominant reactions that can liberate neutrons are 13 C(a, n) l6 0 and 
22 Ne(a, n) 25 Mg. In these reactions, the neutron-rich isotopes, 13 C and 22 Ne 
give up their excess neutrons to heavier nuclei. At this point, we may ask where 
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these excess neutrons came from in the first place. The answer to this interest- 
ing question illustrates an important point about the overall nuclear evolution 
of the universe. 

The abundances that emerge from the Big Bang are roughly 90% by number 
] H and 10% 4 He (e.g. Walker et al 1991). This yields T e = 0.88. On the other 
hand, we may make the observation that and 3 He are the only proton-rich 
(that is, with proton number greater than neutron number) stable isotopes in 
nature. This means that in order for nature to put the nucleons in the universe 
into nuclei with the strongest binding energy per nucleon (iron-group nuclei), 
the Y c of the universe must decrease. 

Most of the decrease in Y c comes from the weak decays in the p-p chains 
and the CNO cycle during hydrogen burning. These interactions drop Y c from 
0.88 to 0.5 in material that has completed hydrogen burning. 4 He itself does 
not have any excess neutrons, but some production of excess neutrons occurs 
in the CNO cycle due to reactions like l2 C(p, y) 13 N(/? + ) 13 C. The net result is 
the conversion of a free proton into an excess neutron, and a drop in Y c . The 
22 Ne production builds up from abundant 14 N produced in the CNO cycle. The 
sequence is 14 N(of, y) l8 F(/? + ) ,8 0(a, y) 22 Ne. Here it is the fact that the only 
stable isotope of flourine is neutron rich that leads to a decrease in Y c . We 
see that the excess neutrons in 13 C and 22 Ne are a consequence of the overall 
drive to decrease Y c in stars. We must keep the goal of the nuclei in sight to 
understand where the excess neutrons come from that drive the s-process. 

The first attempts to understand the details of the s-process led to the classical 
model. The neutron density is always low in the s-process (compared to the 
r-process). If a nucleus is unstable to decay following neutron capture in 
the s-process, it will almost always decay to the first available stable isobar 
before it can capture another neutron. Thus, it generally suffices in s-process 
studies to follow only the abundances as a function of mass number, which 
only change by neutron capture. In this approximation, the rate of change of 
the abundance N A of nuclei with mass number A is 


dN A 

dt 


~n n (av) A N A +n n (av) A -[N A - U 


(9) 


where n n is the neutron number density and {av) A is the thermally averaged 
neutron-capture cross section for the stable isobar of mass number A. We 
can write (cru)^ as a A vj, where Uy is the thermal velocity of neutrons and a A 
is an average cross section, given in terms of i>y. With the definition of the 
neutron exposure 


r = 


J n„v T dt, 


we find 


dN A 


-o a N a + g A -\N A _\. 


( 10 ) 


dr 


( 11 ) 
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Note that the neutron exposure r is a fiuence. It has units of inverse millibarns 
( 1 barn = 10" 24 cm 2 ). Because it is a neutron flux integrated over time, it is an 
appropriate evolutionary parameter for the s-process. If the s-process achieves 
a steady state, then dN A /dx — > 0 and cr A N A -» constant. 

Clayton et al ( 1 96 1 ) were able to show that a single neutron exposure r could 
not reproduce the solar system’s abundance of s-only nuclei. Seeger et al (1965) 
showed that an exponential distribution of exposures, given by 

p(r) = (12) 

TO 

where / is a constant and is the initial abundance of 56 Fe seed, did reproduce 
the solar distribution of s-nuclei. For the distribution of exposures given in 
Equation (12), Clayton & Ward (1974) found that for an exponential average 
of flows in the s-process 

o a n a = fN 56 T 0 Ei 1-1 " 03) 

A' =56 

A fit to the empirical a A N A for s-only nuclei then gives the quantities / and r 0 . 

A complication to the above classical model is the branching that occurs at 
certain isotopes. Here it may be that the decay rate is not considerably greater 
than the neutron-capture rate. In some cases the nucleus may decay before 
neutron capture and in others it may neutron capture before suffering decay. 
The assumptions leading to Equation (9) thus break down. Ward et al (1976) 
developed an analytic treatment of branching in the case of a time-independent 
neutron flux. For time-dependent neutron fluxes, it is necessary in general to 
solve a full network of nuclei numerically (e.g. Howard et al 1986). Since 
the s-process branchings will in general be temperature and neutron density 
dependent, s-nuclei branchings are important diagnostics of the environment in 
which the s-process occurred. We will see this in more detail in Section 4.3. 

4.2 s-Process Sites 

To obtain a good fit of the a N curve to the solar system s-process abundance 
distribution, three distinct exponential distributions of neutron exposures may 
be necessary (Clayton &Rassbach 1967, Clayton & Ward 1974). One exposure, 
with io ~ 0.30 mb" 1 , produces most of the nuclei in the mass range 90 < 
A < 204. This is the main component. Another exposure, with r 0 ^ 0.06 
mb -1 contributes to the A ;$ 90 s-nuclei abundances. This weak component 
is required in order to explain the oN curve around A ~ 90. These two 
components indicate that two separate sites contributed to the abundance of solar 
s-nuclei. Finally, a strong component, with ro ^ 7.0 mb" 1 , may be necessary to 
explain the abundances of the A = 204-209 nuclei. One possible explanation 
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of this component is that the distribution of exposures in the main component is 
not exactly exponential, but rather is higher than exponential at large r. There is 
probably no need for a separate site for the strong component of the s-process. 

The weak s-process component likely comes from He burning in the cores 
of massive stars (£ 15M 0 ) (Truran & Iben 1977, Lamb et al 1977), where 
the temperature is high enough for the 22 Ne(or, «) 25 Mg reaction to produce a 
substantial amount of neutrons. These stars also have strong winds that eject 
this material into the interstellar medium. Recent work has confirmed the 
plausibility of this site (Arnett & Thielemann 1985, Busso & Gallino 1985, 
Prantzos et al 1987, Langer et al 1989, Raiteri et al 1991a, Baraffe et al 1992). 
Uncertainties in the 22 Ne(o?, «) 25 Mg and 22 Ne(a, y) 26 Mg reaction rates prevent 
us from predicting the neutron exposure in these models to high accuracy. 
Recent results on these rates may indicate that the s-process is somewhat more 
robust in this site than previously thought (e.g. Baraffe & El Eid 1994). This 
may complicate the separation of the A ;$ 90 s-nuclei into those coming from 
the weak and main components. 

Some s-processing may also occur in core carbon burning or shell helium 
burning in massive stars. This has been studied by Arcoragi et al (1991) and 
Raiteri et al ( 1 99 1 b). The results indicate that this processing does not contribute 
in a significant way to the weak component. 

The main component of the s-process is likely to occur in the helium-burning 
shell in asymptotic giant branch (AGB) stars (Weigert 1966, Schwarzchild & 
Harm 1967, Ulrich 1973). The structure of such a star is an inert carbon-oxygen 
core, on top of which lies a convective helium-burning shell. On top of this 
helium-burning shell is the hydrogen-rich envelope, which itself is convective. 
The original idea was that the convective helium shell might reach out far 
enough into the hydrogen-rich envelope that protons and 12 C (the result of 
helium burning) could mix and produce I3 C, as discussed in Section 4. 1 . The 
l 3 C would then be the source of neutrons for the s-process. [The current picure 
is that convection does not provide the mixing, but that protons reach down into 
the carbon-rich shell by diffusion or semiconvection (see below).] 

An attractive feature of this model is the fact that the helium burning occurs 
in pulses. Between pulses, hydrogen burns quiescently in a thin shell. Once 
the supply of helium from the hydrogen burning builds up, a helium-burning 
pulse occurs. The energy liberated expands the star and shuts off the hydrogen 
burning. After the pulse has occurred, the star settles down again and begins 
hydrogen-shell burning anew. Pulses last of order tens of years while the 
interpulse periods are of order thousands of years. The significance for the 
s-process is that there is an overlap of mass zones experiencing successive 
helium-buring pulses. Ulrich (1973) was able to show that the mixing and 
burning sequence could naturally give rise to an exponential distribution of 
neutron exposures. Alternating overlap of convection zones can carry the newly 
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produced s-nuclei into the envelope (the so-called “third dredge up”). These 
nuclei would then find their way into the interstellar medium via winds or by 
the ejection of the atmosphere in a planetary nebula phase. 

This nice model for the s-process suffered a setback when Iben showed 
that an entropy barrier prohibited mixing of protons into the helium shell (Iben 
1975a,b, 1976). It was then proposed instead that 22 Ne(a, «) 25 Mg be the source 
(e.g. Iben & Renzini 1983). The helium core grows by accreting the ashes of 
the hydrogen-burning shell. The products of CNO burning are 4 He and I4 N in 
that shell, which combine to give 22 Ne early in helium burning, as discussed in 
Section 4.1. The 22 Ne(a, n) 25 Mg reaction then drives the s-process. The pulse 
and mixing that occurs gives an exponential distribution of neutron exposures. 
This model has some difficulties, however. Basically, the shell flashes in most 
AGB stars are not hot enough to liberate most of the 22 Ne neutrons, and the 
massive ABG stars that are hot enough are too rare. This has led workers to 
consider alternative neutron sources in low-mass AGB stars ( M < 3 M 0 ). 

In low-mass AGB stars, the temperature is too low in the helium-burning shell 
for the 22 Ne(a, n) 25 Mg reaction to be the major source of neutrons. Iben & 
Renzini (1982) argued, however, that, despite the entropy barrier to convection, 
semiconvection or diffusion could cause the mixing of protons with 12 C in 
the interpulse period. This produces pockets of 13 C atop the He zones which 
can liberate neutrons during convective ingestion by the next pulse. Recent 
work indicates that this is a promising site for the s-process (Gallino et al 
1988; Boothroyd & Sackmann 1988a, b,c,d; Hollowell & Iben 1988; Kappeler 
et al 1990). In particular, these models seem to give a good fit to the main 
component of the solar oN curve (e.g. Kappeler et al 1990). We must note 
that these s-process calculations are post-processing calculations, which means 
that the neutron density is a parameterized quantity. Even more serious is the 
lack of a demonstrated occurrence of the needed 13 C-rich pocket, which is 
therefore taken on faith at the present time. It remains to be seen whether the 
good agreement with the solar s-process abundances will hold up when the 
s-process calculations are directly coupled to complete stellar models. Such 
coupled calculations may be available in the not-too-distant future. It will 
also be important to include the effects of energy generation by all the nuclear 
reactions on the stellar structure (Bazan & Lattanzio 1993). 

4.3 Constraints on s-Process Sites 

What constraints can help to evaluate the proposed sites discussed in the pre- 
vious section? s-process branchings are the first important constraints. The 
likelihood that a beta-unstable nucleus in the s-process beta decays depends on 
the rate of beta decay compared to the rate of neutron capture. Evidence for 
branching provides information about these rates. In particular, with knowledge 
of the beta-decay rate from laboratory experiments, the degree of branching 
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constrains the neutron capture rate n n (av). Then knowledge of (cru) from the 
laboratory constrains n n , the neutron number density during the s-process. On 
the other hand, if the beta-decay rate is temperature sensitive (e.g. Takahashi 
& Yokoi 1987), branching data yield constraints on the temperature during the 
s-process. Branching data may also yield constraints on the mass density dur- 
ing the s-process through electron capture rates. Finally, branching data can 
constrain the duration of the neutron pulses (Ward & Newman 1978). If the 
pulse period were much shorter than the lifetime of the branching point isotope, 
there would be no branching. Pulses that were too long in duration would allow 
too much neutron capture. 

What do we find for the s-process in nature? For the main component of 
the s-process, the isotopes l34 Cs, 148 Pm, 151 Sm, 154 Eu, 170 Yb, and 185 W are 
branch-point isotopes with potential as diagnostics of the temperatures and 
neutron densities prevailing during the s-process. Beer et al (1984) used 151 Sm, 
170 Yb, and 185 W to find limits on the neutron number density and temperature. 
Uncertainty in the population of the 137 keV isomeric state in 148 Pm during the 
s-process makes conclusions from this isotope difficult. Uncertainties in cross 
sections and abundances limit the usefulness of l34 Cs and 154 Eu. 

As for the mass density, Yokoi & Takahashi (1983) noticed that l63 Dy could 
beta decay in stars, even though it is stable on Earth. In stars, the 163 Dy atom 
is ionized so that in fact the daughter atom 163 Ho would be at slightly lower 
mass. l63 Ho then could either capture a neutron or electron capture back to 
163 Dy. The electron capture rate depends on the density of electrons, which in 
turn depends on the mass density. Beer et al (1985) were able to constrain the 
mass density in the s-process in this way. 

Finally, Beer & Macklin (1988) studied 151 Sm in order to determine a lower 
limit to the duration of the neutron pulse in the s-process. Studies of 86 Kr may 
give an upper limit to the pulse duration (Beer & Macklin 1989). Unfortunately 
the weak component in this region introduces ambiguities into such an analysis. 

The net results of branching studies in the context of the classical model give 
a temperature for the main component of 2. 8-3. 9 x 10 8 K, a neutron density of 
2.3 — 4.5 x 10 8 cm -3 , a mass density of 2.6-13 x 10 3 gcm~ 3 , and a pulse duration 
of greater than 3 years (Kappeler et al 1989). These numbers agree reasonably 
well with those expected from stellar models. A similar analysis for the weak 
component yields a temperature of 1. 8-3.0 x 10 8 K and a neutron density of 
0.8-1. 9 x 10 8 cm~ 3 (Kappeler et al 1989). 

The relatively high temperatures found in this analysis for the main com- 
ponent suggest that 22 Ne(a, /?) 25 Mg is the neutron source for the s-process. 
Howard et al (1986) studied the s-process nucleosynthesis with this neutron 
source. They obtained poor fits to the solar oN curve when they used param- 
eters derived from stellar models. In particular, the average neutron density 
during the pulses was too high to reproduce the correct branchings. Busso et 
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al (1988) have confirmed these results. On the other hand, the high tempera- 
ture 22 Ne source may simply be the last hot part of the neutron burst that was 
primarily from l3 C at lower temperature (see below). 

Let us consider now the evidence from observations of stars. It was the ob- 
servation of technetium in certain red giant stars (Merrill 1952) that showed that 
stars do indeed synthesize elements and led Cameron ( 1 955) to work out many of 
the details of the s-process. Since all isotopes of Tc are unstable, any Tc present 
in the surface of a star must have been synthesized in the interior of the star by 
the s-process and then dredged up to the surface. Recent observations show that 
red giant stars in the solar neighborhood that do have s-process abundance en- 
hancements in their atmospheres do not show the accompanying enhancements 
of 25 Mg and 26 Mg that one would expect from alpha capture on 22 Ne (e.g. Smith 
& Lambert 1986, McWilliam & Lambert 1988). In addition, observations of 
Rb and % Zr constrain s-process branching at 85 Kr and 95 Zr. Astronomers find 
that the s-process occurring in the interiors of the stars observed must be hap- 
pening at low neutron densities (n n ~ 10 9 cm~ 3 ), not the high neutron densities 
characteristic of the 22 Ne(a, w) 25 Mg reaction (e.g. Lambert 1993). 

From this evidence, it appears that l3 C is more promising as the source of 
s-process neutrons, indicating that low-mass AGB stars are probably the site 
of the s-process. Such stars give a low temperature s-process 1.5 x 10 8 K) 
which would seem to contradict the higher temperatures found from the analysis 
of the s-process branchings in the classical model ( T = 2. 8-3. 9 x 10 8 K) 
discussed above. In the low-mass AGB star s-process calculations that do show 
good agreement with solar abundances (e.g. Gallino et al 1988, Kiippeler et 
al 1990), there are two bursts of neutrons per pulse: a strong burst due to the 
13 C(a, n) l6 C reaction at T ~ 1.5 x 10 8 K, and a second, weaker one, due to 
the 22 Ne(ar, «) 25 Mg reaction. This weaker burst occurs when the helium shell 
contracts following the first burst and heats to a temperature of T ~ 3x 10 8 K. It 
resets the branch point thermometers to this higher temperature, in agreement 
with the analysis from the classical model. 

More evidence for 13 C as the dominant source for neutrons in the s-process 
comes from studies of galactic abundance evolution. Mathews et al ( 1 992, 1 993) 
studied the evolution of the Ba/Fe ratio in our Galaxy. Ba is predominantly an 
s-process element and hence must be secondary (i.e. made from initial Fe). 
Mathews et al found that only an s-process behaving as a primary process fit 
well the observations of Ba abundances in the atmospheres of old stars. 1 3 C is a 
primary neutron source, as discussed in Section 3.3. 22 Ne is secondary because 
it must be built up from pre-existing CNO nuclei. The Fe seeds are of course 
secondary. Clayton (1988a) described how the secondary s-process with the 
l3 C neutron source is able to mimic primary nucleosynthesis. The idea here is 
that while the galactic abundance of Fe seed for the s-process grows with time, 
so does the abundance of s-process neutron poisons. 
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A final point of increasing relevance is the new information from pre-solar 
SiC grains found in the Murray and Murchison meteorites. These grains are 
carriers of isotopic anomalies in s-process isotopes (Srinivasan & Anders 1978, 
Tang & Anders 1988). In addition, these grains are anomalous in their Si and 
C (Zinner et al 1987, Anders & Zinner 1993). It appears that these grains have 
condensed in carbon-star atmospheres, which are s-process enriched and have 
variable 13 C-rich compositions (Lambert et al 1986). As surviving stardust, the 
grains are almost pure endmembers in the “cosmic chemical memory” theory for 
interpreting isotopic anomalies in solar system samples (Clayton 1978, 1982). 
From studies of trace s-isotopes in these grains (Ott & Begemann 1 990 a,b; Zin- 
ner et al 1991; Richter et al 1 992; see Anders & Zinner 1 993 for a review) the 
abundance ratios only fit if the grains come from low-mass AGB stars (Gallino 
et al 1990). A vexing problem with this idea, however, is that such stars cannot 
explain the anomalous Si isotopes (a major constituent of the grains). One sug- 
gested answer is higher mass AGB stars, in which burning of Mg isotopes in late 
pulses resets the ratio of 29 Si/ 30 Si(Brown & Clayton 1992). Galactic abundance 
evolution of Si isotopes may also hold the key (Clayton 1988b, Gallino et al 
1994). Alternatively, some other site may be responsible for these grains (e.g. 
Arnould & Howard 1993). These tiny, sturdy grains have traveled from afar car- 
rying important messages about the s-process which have yet to be deciphered. 

In summary, low-mass AGB stars are at present the most promising site for 
the main component of the s-process. Confirmation of this site will require 
continued interplay of nuclear physics, meteoritics, stellar evolution and struc- 
ture theory, nucleosynthesis theory, galactic abundance evolution theory, and 
stellar astronomy. Many people will be busy for quite some time to come! 

5. THE p-PROCESS 


. . . and the elements shall melt with fervent heat . . . 


II Peter 3: 10 

We turn finally to the p-nuclei. These are the 35 nuclei bypassed by the r- 
and s-processes. As we see from Figure 1, except for the light p-nuclei ( 92 Mo, 
94 Mo, 96 Ru, 98 Ru), the abundances of p-nuclei are considerably less than those 
of their r- and s-nuclei counterparts. Furthermore, the p-process abundance 
distribution shows interesting structure with peaks at 92 Mo and 144 Sm. These 
are important clues for determining where the p-process occurs. 

It is probably wrong to think that the p-process occurs in a single site. We 
can imagine many astrophysical settings where conditions are right to modify 
a pre-existing supply of r- and s-nuclei to form p-nuclei. The relevant question 
is really what site contributes the bulk of the p-nuclei. For more details on the 
p-process, the reader should consult the excellent review by Lambert (1992). 
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5 . 1 General Considerations 

It is impossible to produce p-nuclei by neutron capture. How then can nature 
make these nuclei? The first possibility that suggests itself is proton capture. 
It may be that in the course of the evolution of some system striving to reach 
NSE, protons are liberated which can capture on pre-existing seed nuclei to 
make p-nuclei. Alternatively it may be that in a freezeout from proton-rich 
NSE, free protons could capture on seed iron-group nuclei. reactions could 
allow further capture of protons to higher mass. 

We can explore the conditions required for such p-processes with the help 
of Figures 9 and 10. Figure 9 shows the timescale for capture of a proton [a 
(p,y) reaction] by the most proton-rich beta-stable isotope of each element 
at the fixed temperature of Tg = 1 for different values of the proton mass 
density pY p . The timescales for a proton-disintegration (y, p) reaction or for 
a neutron-disintegration (y, n ) reaction for these isotopes are all greater than 
10 10 seconds. Note that for a site with a proton mass density pY p of 1 g cm -3 , it 
would take 92 Mo (Z = 42) about 10 4 s to capture a proton. Capture of protons 
on higher-charge isotopes would take even longer. The timescale for proton 
capture decreases if the setting has a higher density of protons available. For 
example, for pY p = 10 3 gcm~ 3 , the timescale for capture of protons on 92 Mo 



Figure 9 Timescales for proton capture on the most proton-rich isotope of each element at the 
fixed temperature of T 9 — 1 . The curves are for mass densities in protons of pY p = 1 g cm -3 ( solid 
curve), 10 3 gcm -3 (long dashed-dotted curve), and 10 6 gcm -3 (short dashed-dotted curve). The 
rates are computed from expressions in Woosley et al (1975). 
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Figure 10 Same as Figure 9 but for T 9 = 3. The shorl-dashed jagged curve shows the timescale 
for ( y , n) reactions on these nuclei. The (y, n) rates are computed from neutron-capture cross 
sections in Woosley &. Hoffman or Cowan et al (1991 ) and neutron -separation energies are derived 
from Moller& Nix (1988). 

would be about 10s while for pY p = 10 6 gcm~ 3 it would be about 10 _2 s. If 
an astrophysical site could maintain a mass density in protons of 10 6 g cirT 3 for 
1 0 5 s at Tg = 1 , the proton-rich isotopes of all elements up to platinum (Z = 78) 
could capture a proton. The question for finding the p-process site is whether 
such conditions are possible. It is unlikely. 

The proton-capture rates also increase with increasing temperature because 
the reactants have a higher relative kinetic energy compared to the Coulomb 
barrier than at lower temperature. This leads us to ask what happens to the 
timescales if we increase the temperature. We see the effect in Figure 10 for 
To = 3. The timescale for proton capture does indeed decrease, but so does 
the timescale for a (y, n ) reaction. For pY p = 1 gem -3 , the timescale for a 
(y, n) reaction is less than that for a (/?, y) reaction at Z & 40. This means 
that it is more likely under these conditions for a proton-rich nucleus to suffer 
a (y, n) reaction than to capture a proton. Of course this makes sense from our 
discussion in Section 2. If the system is evolving towards NSE, nuclei more 
massive than the nucleus with the highest binding energy per nucleon will tend 
to disintegrate nucleons to increase the number of macroscopic states available 
to the system. 

For higher proton mass densities, the (y, n) reactions do not dominate the 
proton captures until higher nuclear charge: Z % 40-50 for pY p = 10 3 gcm 3 
and Z as 70 for pY p — 10 6 gcm -3 . One might imagine that the system could 
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then produce p-nuclei at least up to ytterbium (Z = 70) under such high tem- 
perature and high density conditions. Such conditions are extremely difficult 
for nature to achieve, however. The dilemma for making p-nuclei is clear. If 
nature is to make these nuclei by proton capture at low temperature where the 
flow will not be impeded by disintegrations, a large supply of protons must be 
available for a long time. If the proton capture process is to occur at higher 
temperature where the capture timescales are shorter, disintegration reactions 
will dominate the flow and prevent capture to higher mass. 

The escape from this dilemma is the realization that proton captures need 
not make most of the p-nuclei. The various disintegration reactions do the job. 
In particular, if pre-existing r- and s-nuclei are exposed to high temperature, 
nuclear reactions will occur and tend to drive the abundances toward NSE. 
The first reactions to occur are the (y, n) reations which produce quite proton- 
rich nuclei. Once the nuclei become sufficiently proton rich, they then begin a 
(y, p) and (y, a) cascade. In this way the nuclei “melt” towards iron. If the high 
temperature drops off quickly enough, the system does not reach NSE and the 
melting will be incomplete, leaving an abundance of proton-rich heavy nuclei — 
the p-nuclei. Where the disintegration flow crosses the N = 50 and N = 82 
closed neutron shells, the disintegration timescales become large because of the 
particularly strong binding energies. This is apparent in Figure 10 as the peaks 
in the (y, n ) timescales at the closed shell nuclei 92 Mo (Z = 42) and ,44 Sm 
(Z = 62). Because of the long disintegration timescales, abundances build 
up at these nuclei. In this way we can explain the peaks in the solar system’s 
p-process abundance distribution. As a final note, the extremely short (y, n) 
timescale for 180 Ta (Z = 73 in Figure 10) explains why this fragile species is 
the rarest stable isotope in nature. 

5.2 p-Process Sites 

The early papers on the p-process considered the proton capture mechanism 
(Burbidge et al 1957, Ito 1961, Macklin 1970, Truran & Cameron 1972, Au- 
douze & Truran 1975). The site was imagined to be the hydrogen-rich enve- 
lope in massive stars undergoing a supernova explosion. The supernova shock 
passing through this region would heat up the material and proton capture 
reactions would produce the p-nuclei. However, the densities, temperatures, 
and timescales required are unrealistic for the hydrogen-rich envelope (see, for 
example, the discussion in Woosley & Howard 1978). 

Arnould (1976) computed the p-process in the hydrostatic oxygen burning 
phase in stars. The timescales are longer in this site than in the supernova 
site and would allow for more proton capture. In this site, temperatures were 
high enough for disintegrations [especially (y , n ) reactions] to be important. A 
major challenge for this model is to eject the new p-nuclei without significantly 
modifying their abundances during the subsequent supernova explosion. 
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Woosley & Howard (1978) computed the p-process in the O/Ne shell in type 
II, that is, core-collapse, supernovae. The supernova shock heats up this shell 
and causes the partial melting of the nuclei. In this model, only disintegra- 
tions are important, hence the alternative name “gamma-process” (see Rayet 
et al 1990, Prantzos et al 1990, Rayet et al 1992 for extensions of this model). 
This model successfully reproduced most of the p-nuclei in their solar system 
proportions, although it seriously underproduced the light p-nuclei. As in the 
hydrostatic oxygen burning model, it was necessary to superimpose several 
abundance distributions to get a realistic distribution of p-nuclei. Type II su- 
pernovae should naturally give a distribution of conditions depending on the 
layers of the proto-SN considered that would naturally give a distribution of 
abundances. The inner regions of the O/Ne shell will achieve the highest tem- 
peratures and thus get closest to NSE. These regions make the lighter p-nuclei. 
Outer regions produce the heavier p-nuclei because the “melting” is less com- 
plete. Prantzos et al (1990) computed the p-process abundance distribution for 
a specific model of supernova 1987 A and found that a distribution of condi- 
tions naturally arose and gave a solar system p-process abundance distribution, 
except for underproduction of light p-nuclei. 

The problems with the underproduction of the light p-nuclei in the gamma- 
process led Howard et al (1991) to consider the production of p-nuclei in the 
outermost layers of a carbon-oxygen white dwarf star suffering a type la su- 
pernova explosion. In this model, s-processing prior to the explosion built up 
the abundances of A & 90 nuclei. The high density during the explosion then 
allowed proton capture reactions to produce many of the light p-nuclei while 
the normal gamma-process made the heavier p-nuclei. Later calculations using 
the realistic type la models of Khoklov (1990) have not been as successful in 
producing the light p-nuclei (Howard & Meyer 1992). More studies of this 
promising site are required. 

Some p-nuclei may also be produced in spallation reactions. Most notably, 
neutrinos may spall neutrons from heavy nuclei during type II supernovae to 
make p-nuclei. Such a process is only likely to produce significant amounts of 
the rarest p-nuclei such as 138 La and 180 Ta (Woosley et al 1990). 

5.3 Some Constraints on p-Process Models 

As mentioned above, the p-process probably occurs in several places in nature. 
Any astrophysical setting in which high temperatures but sufficiently short 
timescales lead to incomplete melting of heavy nuclei can produce p-nuclei. 
What we really seek is the site that produces the bulk of the p-nuclei. We can 
attempt to analyze this requirement by means of overproduction factors. 

The overproduction factor of an isotope in the product material of some 
nucleosynthetic process is the ratio of its mass fraction in that material to its 
mass fraction in the solar system. In order for the process to be responsible 
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for the production of the bulk of the solar system’s supply of a given isotope, 
that isotope must have the largest overproduction factor in the product ma- 
terial. If the process is to be responsible for the bulk of the solar system’s 
supply of two or more isotopes, those isotopes must all have comparably large 
overproduction factors. 

Prantzos et al (1990) found an average overproduction factor of 0.96 for 
p-nuclei in 15M 0 of ejecta in their model for SN 1987A. The overproduction 
factor for 16 0 in the same model was 1 1 .5 (Thielemann et al 1990). Since type 
II supernovae made most of the solar system’s 16 0, this would indicate that the 
type II supernova site could not be responsible tor the production of the bulk ot 
the solar system p-process elements. Prantzos et al argue, however, that there are 
extenuating factors. First, the SN 1987 A model used Large Magallenic Cloud 
metallicities, which are lower than those in our Galaxy. Milky Way metallicities 
could give enhanced seed abundances and thus higher p-process overproduction 
factors. Also, it may be that p-nuclei production relative to oxygen could be 
higher in higher mass stars. Surveys over a range of star masses, such as that 
of Arnould et al (1992), will be important for understanding the contribution 
of type II supernovae to the solar system p-process abundances. 

As for the type la model, Howard et al (1991) noted that besides possibly 
producing p-nuclei, these supernovae made most of the solar system’s 56 Fe. 
The requirements that type la supernova models make 56 Fe and p-nuclei in 
solar proportions and that they produce 0.5-1 .0 M Q of 56 Fe, and the fact that the 
typical p-process overproduction factors in this model are ~ 10 , lead to the 
conclusion that the zones that produce p-nuclei comprise 0.04-0.08M Q of the 
white dwarf. This is in good agreement with the models (e.g. Khoklov 1990). 
It appears that type la supernovae are capable of producing the bulk of the solar 
system’s p-nuclei. 

Another important constraint on the p-process is the presence of live ’Sm 
in the early solar system (Lugmair et al 1983, Prinzhofer et al 1989). From 
general galactic abundance evolution arguments, Prinzhoter et al inferred from 
their measurements that the production ratio ot l46 Sm/ 144 Sm should be between 
0.07 and 0.5. This ratio causes problems for the gamma-process for which the 
production ratio for these isotopes is typically —0.02 tor the type II model 
(Woosley & Howard 1978) and -0.05 for the type la model (Howard et al 
1991). The inferred production ratio is in fairly good agreement with that 
found from proton-capture models (e.g. Audouze & Truran 1975), which led 
Prinzhofer et al to tavor such models tor production of p-process elements. 
As we have seen, however, such models are not astrophysically realistic, so 
the measurements present a challenge for gamma-process models. Woosley 
& Howard (1990) found, however, that if the branching ratios for (y, n) and 
(y (x ) reactions on * 4 ^Gd are varied within experimental uncertainties, the 
,46 Sm/ 144 Sm production ratio could be increased dramatically. It is clearly 
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a worthy goal for experimental nuclear physicists to attempt to make accurate 
measurements of the disintegration rates of 1 48 Gd. Another important effect that 
may alleviate the 146 Sm/ 144 Sm production problem is galactic infall (Clayton 
et al 1993). 

In summary, plausible models exist that do a good job of satisfying most 
ot the constraints on the p-process. The most vexing puzzle that remains is 
the underproduction of the light p-nuclei. Perhaps, with more work, we will 
see that type la supernovae can produce these nuclei at the correct levels. On 
the other hand, it may be that we will need to turn our attention to some other 
site, such as the ^-process in the mass cut in type II supernovae (Woosley & 
Hoffman 1992) or proton-capture reactions in Thorne-Zytkow objects (Cannon 
et al 1992), to explain the origin of these isotopes. 


6. CONCLUSION 

The terms r-, s-, and p-processes have been successful in clarifying our picture 
of heavy-element formation. In this paper, we have seen that this is because 
these different processes reflect the responses a nucleosynthetic system can have 
to being out of equilibrium. In this way, our understanding of the mechanisms 
of heavy-element synthesis is firm. 

Future work will continue to focus on using all clues at our disposal to resolve 
the problems surrounding the astrophysical sites for these processes. Important 
questions we must seek answers to in the near future include: 1. Do nascent 
neutron star winds have high enough entropies to drive a full r-process?, 2. 
Do I3 C pockets really form in low-mass AGB stars, and do they give the right 
neutron exposures for the s-process?, and 3. Can type la supernovae make the 
p-nucleus 92 Mo? That we can ask such specific questions indicates that we 
have good ideas about the sites for the r-, s-, and p-processes. That we do not 
have answers shows that we have much work to do. 
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